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Abstract  of  Dissertation  Presented  to  the  Graduate  Council  in 
Partial  Fulfillment  of  the  Requirements  for  the  Degree  of 
Doctor  of  Philosophy 

GENERAL  NONLINEAR  THEORY  OF  SANDWICH  SHELLS 

by 
Ju-Chin  Huang 
December,  1965 
Chairman:  Dr.  I.  K.  Ebcioglu 
Major  Department:  Engineering  Science  and  Mechanics 

A  general  nonlinear  theory  of  sandwich  shells  to  the  full 
extent  of  the  nonlinear  strain-displacement  relations 

^  '^.j  =  JC  Vtij-^Vjic-^  V^WIj  ) 

,   has  been  obtained  by  means  of  the  modified  He Hinge r-Reissner 
variational  principle  of  three-dimensional  elasticity.   The 
fundamental  equations  are  in  tensor  notation  and  in  terms  of  the 
undeformed  state.   By  this  same  technique  and  by  making  certain 
simplifying  assumptions  in  the  strain-displacement  relations,  par- 
tially nonlinear  theories  for  sandwich  shells  are  also  obtained. 
These  approximations  are  based  on  the  assumption  of  small  strain 

I  but  large  deflections,  thin  facings,  and  soft  cores,  respectively. 

These  intermediate  theories  are  more  suitable  for  application. 

The  intermediate  and  the  linearized  form  of  the  various  sets 
of  equations  derived  herein  coincide  with  known  theories. 
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CHAPTER  I 


INTRODUCTION 


A  sandwich  structure  is  formed  by  two  thin  facings  of  a 
strong  material  between  which  a  thick  layer  of  very  light-weight 
and  comparatively  weak  material  is  sandwiched.   The  advantage  of 
this  kind  on  construction  is  the  large  moment  of  inertia  of  the 
section  provided  by  spacing  far  apart  the  load-carrying  facings. 
Accordingly,  the  applications  of  sandwich  structures  in  various 
areas,  especially  by  the  aerospace  industry,  have  increased  in 
recent  years.  As  a  result,  more  research  work  concerning  this 
kind  of  structure  has  become  desirable.   A  series  of  extensive 
bibliographies  (1),  (2)*  and  comprehensive  reviews  on  the  analysis 
of  sandwich  structures  by  Habip  (3),  (4)  have  appeared  recently. 

In  the  early  work  of  Reissner  (5),  (6),  (7)  and  Wang  (8), 
the  simplest  model  consisting  of  two  facings  acting  as  membranes 
and  a  core  resisting  transverse  shear  and  normal  stresses  has 
been  employed  for  deriving  the  governing  equations  of  sandwich 
structures.   In  recent  years,  it  appears  that  changes  in  technol- 
ogy and  a  concern  for  the  optimization  of  structural  elements 


*Under lined  numbers  in  parentheses  refer  to  the  Bibliography 
at  the  end  of  this  dissertation. 


» 


> 


subjected  to  thermal  as  well  as  mechanical  loads  have  brought  about 
several  studies  of  a  new  type  of  sandwich  construction  with  strong 
cores.   The  theory  then  takes  the  flexural  rigidity  as  well  as 
transverse  normal  deformation  of  the  core  into  account  while 
including,  as  usual,  the  flexural  rigidities  of  the  upper  and  lower 
facings  about  their  own  middle  surfaces.   Grigoliuk  and  Chulkov  (9) 
have  presented  a  paper  on  this  subject  for  the  case  of  small  de- 
flection theory  of  sandwich  shells.   They  consider  the  core  as  a 
three-dimensional  body  and  assume  that  the  displacements  can  be 
expressed  approximately  as  a  linear  function  of  the  transverse 
coordinate.  Wu  (10)  generalized  this  to  a  large-deflection  theory 
of  orthotropic  sandwich  shallow  shells.   In  the  latest  works  of 
Ebcioglu  (11) ,  (12) ,  the  nonlinear  field  equations  for  the  sandwich 
plates  and  shells  have  been  obtained  by  means  of  the  Hamilton 
principle.   For  plates,  the  nonlinear  strain-displacement  relations 
are  used  to  the  full  extent,  but  partially  nonlinear  strain-displace- 
ment relations  are  used  to  the  full  extent,  but  partially  nonlinear 
strain-displacement  relations  only  are  employed  for  shells.   The 
present  dissertation  is  inspired  by  Habip's  works  (13),  (14) .   The 
fundamental  equations  of  the  theory  of  plates  and  shells  have  been 
obtained  by  Habip  from  the  three-dimensional  theory  of  nonlinear 
elasticity  by  integration  across  the  thickness  of  the  undeformed 
plate  and  by  means  of  a  modified  version  of  the  Hellinger-Reissner 
variational  theorem  of  three-dimensional  continuum  djmamics, 
respectively,  for  the  case  when  the  "shifted"  components  of  dis- 
placement can  be  assumed  to  vary  linearly  through  the  thickness  of 
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the  shell. 

The  present  study  is  an  attempt  at  obtaining  a  general  non- 
linear theory  of  sandwich  shells  to  the  full  extent  of  the  general 
nonlinear  strain-displacement  relations  (39)  in  tensor  notation  and 
in  terms  of  the  undeformed  state  by  means  of  the  modified  Hellinger- 
Reissner  variational  principle  of  three-dimensional  elasticity. 
The  method  employed  here  is  similar  to  that  used  by  Habip  but  the 
order  of  variation  and  "shifting"  is  slightly  different.   In  this 
work,  we  introduce  the  geodesic  normal  coordinate  system  into  the 
variational  equation  before  we  carry  out  the  variation.   In  the 
general  theory,  the  results  are  identical  regardless  of  the  order 
of  variation  and  "shifting"  but  in  partially  nonlinear  theories  the 
results  come  out  quite  different.  We  shall  discuss  this  in  more 
detail  in  Chapter  V. 

In  Chapter  III,  the  "exact"  fundamental  equations,  in  the 
sense  of  using  the  complete  general  nonlinear  strain-displacements 
relations,  are  given.   No  assumptions  about  the  state  of  deformation 
have  been  made  except  that  the  displacements  vary  linearly  through 
the  thickness.   But  these  "exact"  equations  are,  for  practical 
purposes,  too  complicated;   so  we  introduce  some  simplifications 
in  Chapter  IV.   These  approximations  are  based  on  the  assumption 
of  small  strain  but  large  deflections,  thin  facings,  and  a  soft 
core,  respectively.   By  making  these  simplifications  we  arrive  at 
several  approximate  theories  suitable  for  applications.  When  these 
theories  are  compared  to  some  known  results,  they  do  agree  exactly. 


A^;'r"% 


The  present  work,  as  in  (11) ,  (12) ,  takes  into  account  the 
effects  of  transverse  shear  and  normal  stresses  as  well  as  rotatory 
inertia,  with  different  material  densities  and  material  constants 
in  each  layer  of  the  sandwich  shell.   In  addition,  each  of  the 

I  three  layers  is  of  different  thickness,  and  no  a  priori  limitations 

are  imposed  upon  the  displacement  functions  until  Chapter  IV,   It 

,  •  r        is  also  assumed  that  the  facings  and  the  core  are  anisotropic, 
having  elastic  symmetry  with  respect  to  the  middle  plane  of  the 

'.  •  layers.  The  effect  of  steady  thermal  gradients  is  also  included 

in  the  stress-strain  relations. 
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CHAPTER  II 
PRELIMINARIES 


In  order  to  formulate  a  geometrically  nonlinear  theory  of 
sandwich  shells,  some  fundamental  concepts  from  tensor  analysis, 
surface  geometry,  and  elasticity  will  be  used  repeatedly  in  this 
study  and  are  reproduced  here  for  convenient  reference.  More 
detailed  treatment  of  these  subjects  may  be  obtained  in  references 
(15),  (L6),  and  (17). 

2.1  Outline  of  Tensor  Analysis 

^'      Convected  coordinates  and  base  vectors.   Let  £ '  be  a 
set  of  fixed  rectangular  cartesian  coordinates.   The   €,•   are 
unit  base  vectors,  and  the  position  of  a  point  in  space  with 
coordinates   Z  *•   can  be  defined  by  the  position  vector  f 


where 


The  range  of  Latin  indices  is  1,  2,  3.  Repeated  indices  are  to  be 
summed  over  their  range. 

The  differential  of  the  position  vector  in  (1)  is 

dr=lL,dz'=jiz'e,       .  (2) 
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The  length  of  a  line  element  is  defined  as 


(3) 


where  hcj      denotes  the  Kronecker  symbol. 

Let  us  now  introduce  a  general  convected  coordinate  system 
©     ,  defined  by  the  coordinate  transformation 


9"=  e'CH',  z^  z'; 


(4) 


or 


r=  z'ce\  e\  e'  ) 


(5) 


provided  that  the  Jacobian  of  (4)  does  not  vanish. 

If   r    is  regarded  as  a  function  of   9    ,  and  using 
a  comma  followed  by  a  subscript  i.   to  denote  partial  differentiation 
With  respect  to  9    , 

where 


9i  =  r,o    = 


d9' 


(7) 


are  the  base  vectors  of  the  convected  coordinates. 

By  the  chain  rule  of  partial  differentiation,  the  differential 
o  "     can  be  expressed  as 


> 

^ 


> 


Substitution  of  this  expression  into  Equation  (6)  yields 
which,  by  comparison  with  (2),  leads  to 

Sj  If^   =    f t  •  (w) 

From  the  definition  of  the  law  of  covariant  transformation  we  know 
that    ^  c        are  covariant  base  vectors  for  the   Q  coordinate 

system. 

By  using  (3),  the  length  of  a  line  element  is 

J^"=cif-clf  =  ^^.^.de'<ie'=  ^.jde'de'  ^  (ii) 

where 

is  called  the  covariant  metric  tensor  for  the    0     coordinate 
system.   From  Equations  (3),  (11),  and  with  the  help  of  Equation 
(8),    . 

^''        2z'    ^z^    ^•""   • 

B.   Reciprocal  base  vectors.   There  is  no  distinction  between 
covariant  and  contravariant  components  of  the  base  vectors   Gi 
in  a  Cartesian  coordinate  system.  Thus  we  can  write  6*  =  fi^    , 


I 
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and  define  contravariant  base  vectors  in  the   0    coordinate 
system  by  the  law  of  contravariant  transformation  (17) ,  namely. 

In  view  of  Equation  (10b), 

3    •  3j  =  ^j      •  .  (15) 

so  that,  in  general,  the  covariant  and  contravariant  base  vectors 
in  the   0    coordinate  system  are  orthogonal  and  reciprocal  to 
each  other.  By  analogy  with  Equation  (12),  the  contravariant 
metric  tensor  is  defined  by 

9''=  3^-  r  .  (16) 

,  The  following  useful  formulae  are  derived  from  Equations 
(14),    (10),   and    (15).  ' 

3'=^3''3j  ■  .  (17b) 

^        9jA    =     ^  J         .  .  (17c) 

Thus,  the  metric  tensor  can  be  used  to  raise  or  lower  indices. 
The  volume  element  for  the  convected  coordinate  system  0         is 
given  by 

Jr  =JYcie'66^de^  ,  (is) 


where 


3  = 


(19) 


> 


f^\ 


C.  Derivative  of  a  vector.   The  covariant  base  vectors 
^ i,  are  expressed  in  terms  of  the  position  vector   f    by 

Equation  (7) .   Since  f        is  assumed  to  be  a  continuous  function 


of   0   ,  it  follows  that 


g^>j  =  3^>i 


Differentiating  Equation  (10a),  and  using  (10b),  we  obtain 


(20) 


9m  j  =^'"-:r7I 


-^-Z""  ^^1  o-r 


Introducing  the  definitions  of  Christoffel  symbols  of  the  first 
and  second  kind,  respectively. 


I  oj  r      Omn 


TZ*"    "^H* 


■^e"  je-'  Be*" 


and 


(21) 


(22) 


the  derivative  of  base  vectors  (20)  becomes 

Now  for  any  vector  \/   ,  which  is  expressed  in  terms  of 
its  components  by 

v  =  Vir  =  v'fc  ,  (23) 


^' 


> 


10 
the  derivative  can  be  written  as 

v,j  =  (^^!i^.j  =  v\j9i^v'!i,j     .  (24) 

By  using  Equation  (22) , 

V.j  =  V-lj  9«      ,  (25)  , 

where 

vHj  =  v'.j  +  r".  v'  (26) 

is  the  covariant  derivative  of  the  vector  V 
Similarly, 

V.j  =  /-Ij  f "  (27) 

where 

Vtlo=  ^-j  -^rjV.  (28) 

is  the  covariant  derivative  of  the  vector  ^i       . 

2.2  Concepts  from  Three-dimensional  Theory  of  Elasticity 

A.   Strain  tensor.   Let  the  undeformed  body  be  described 
by  a  general  right-handed  convected  coordinate  system.    © 
If  the  position  vector  of  a  point  in  the  undeformed  body  is 
denoted  by 

f  =  Tce\  d\  e" )      ,  (29) 

then  when  this  body  is  subjected  to  load,  it  will  deform  into 
a  new  configuration 

R=R(e',  G\e')  (30) 

which  is  the  position  vector  of  a  point  in  the  deformed  body. 


11 

The  0      are  material  coordinates  which  are  associated  with 
corresponding  points  in  the  deformed  and  undeformed  body. 

In  the  undeformed  body,  the  length  of  a  line  element  is 
given  by 

J^^=   5o  de^c^e^'  .  (31) 

In  the  deformed  body,  the  length  of  the  same  line  element  becomes 

d$'=  Gi,j<i9'de'  .  (32) 

where  CriJ    is  the  metric  tensor  in  the  deformed  state.   The 
difference  in  lengths  of  these  line  elements  is  a  measure  of  the 
deformation  experienced  by  the  body  as  it  moves  from  the  undeformed 
to  the  deformed  position.  This  deformation  is  described  by  the  strain 
tensor  defined  by 

dS^-clA^=    2r.jJl9'd9^  .  (33) 

from  which  it  follows  that 

This  shows  that  oi:    is  a  symmetric  tensor. 

The  strain  tensor  can  be  expressed  in  terms  of  the  displace- 
ment  vector  V     by  noting  that 

R  =  r  +■  \7  .  ■   (35) 

By  using  the  formulae  given  in  the  previous  section,  the  base 
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vectors  and  metric  tensor  for  the  convected  coordinates  0   may 
be  defined  in  both  the  undeformed  and  the  deformed  body.   Thus, 


(36) 


In  view  of  Equations  (27)  and  (35) ,  the  base  vectors  Cf^    and 
metric  tensor  ct^j  become 

5,=(sr+ v"/. )  ?^  ,  (37) 

c^.,  =  (ST.rucS%v"g5_         .  (38) 

Then,  substituting  these  into  Equation  (34), 

rij=j(Vjj^V;/.  +  V^V^rij)         .  (39) 

B.   Stress  tensor.  The  state  of  stress  is  here  defined  in 
terms  of  the  three-dimensional  stress  vector,  ^X       >  per  unit  area 
of  the  undeformed  body  and  associated  with  a  surface  in  the  deformed 
body  whose  unit  normal  in  its  undeformed  position  is  ,n   ,  namely 
(18), 

^-^''>^^^^  ,  (40) 

where 


> 
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By  using  Equation  (37),  expression  (40)  becomes 

.t  =  s'^StVlJ.n,  g,  =  t^\^.  £.     .  (41) 

and  T   are  contravariant  stress  tensors  measured  per 
unit  area  of  the  undeformed  body  and  referred  to  base  vectors  in 
the  deformed  and  undeformed  body,  respectively. 

C.   Stress-strain  relations.   The  conventional  Duhamel-Neumann 
stress-strain  relations  for  an  isotropic  material  can  be  generalized 
for  an  anisotropic  material  as  follows  (13) . 

^M=£g^nS'""-H-cjT  .  (42) 

The  inverse  of  (42)  is 

S''=     C''''(     Yrs-^r.-r)  .  (43) 

The  coefficients  P..     /-'•J'"    aj-g  the  elastic  coefficients  of 
the  medium.   They  depend  on  the  metric  tensors  and  physical 
properties  of  the  undeformed  body.  The  oC;j   are  the  coefficients 
of  thermal  expansion.  All  of  these  coefficients  are  also  functions 
of  0      through  the  steady  temperature  field  T  and  satisfy  the 
following  symmetry  conditions 

P     =■?     =P,.=P.- 

^-Cjmn      '-Jimn     '-i.jnm  '-/nnij       , 
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For  a  medium  having  elastic  symmetry  with  respect  to  the  surface 
9   =■     const.,  all  coefficients  containing  the  index  3  either  once 
or  three  times  vanish  (15) ,  i.e., 

C^^^=   C'"  =  O       .  (44) 

Equations  (43)  then  reduce  to  (14) 


(45a) 


5"'=  2  C'''/'^CT/,3-<>v».T)   , 


7»»  '  -^  >  (45b) 

5   =C      C"J(Sv-0^/>.T)+  C     C^3>-C^,3T;   .       (45c) 


''■'.  2.3  Surface  Geometry 

A  shell  is  a  body  occupying  the  space  between  two  surfaces 
(called  faces)  a  small  distance  apart.   The  coordinate  system  is  so 
chosen  that  the  surface  defined  by  0  =  o  lies  midway  between  the 
faces.   It  is  called  the  middle  surface.  When  B       is  measured 
along  a  line  perpendicular  to  the  middle  surface,  the  coordinate 
;  "  :        system  is  called  normal  (17) .   In  a  normal  coordinate  system  the 
position  vector  f  for  the  undeformed  body  becomes 

The  vector  T   locates  points  on  the  middle  surface,  and   O 
and  y   are  general  curvilinear  coordinates  on  this  surface.   The 
vector  a J  is  the  unit  normal  to  the  middle  surface,  here  directed 
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> 


outward  on  a  surface  of  constant  positive  Gaussian  curvature. 
The  covariant  base  vectors  of  the  middle  surface  are 


oc 


(47) 


Here  and  in  what  follows,  we  use  the  convention  that  Greek  indices 
take  on  the  values  1,  2.  The  metric  tensor  of  the  middle  surface  is 


The  contravariant  components  are  defined  by 


(48) 


a''^  CL^^  =  I 


(49) 


The  contravariant  base  vectors  are  defined  by 


r> 


The  area  of  an  element  on  the  middle  surface  is 


(50) 


J/\  =  /^  <yd'  ^e" 


(51) 


where 


CL  = 


a. 


=  a„  dx,.-Ca,i) 


(52) 


The  coefficients  of  the  second  quadratic  form  are  defined  by 


La ^w  ■  ^J,/5  —  ^3  •  ^~, 


f 


(53) 


The  mixed  components  of  the  second  fundamental  form  are 


\)\  -     (^"^  ^> 


(54) 
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and  the  contravariant  components  are 

\>'''=  a"nf  •  (55) 

The  Christoffel  symbols  for  the  middle  surface  can  be  obtained 
by  evaluating  Equation  (22)  at  0  =  ^  .  We  have 

a.,^  =  r'JJ^  ^r  +  r?^  ^,        ..  ,  (56) 

where  an  asterisk  denotes  quantities  on  the  middle  surface. 
Equation  (56)  may  be  rewritten  as 

docWp,  =    \     '^/h    ^i  .  (57) 

Here  the  double  vertical  line  stands  for  covariant  differentiation 
based  on  the  metric  tensor  (^o/a  • 

From  Equations  (53)  and  (21),  when  evaluated  at  &  "  o    ^ 


and  hence 


(58) 


Ci^llfl,=    \>u/i   a,  .  (59) 


\\p>~    o,^^ 


Multiplying  Equation  (53)  by  dX,   ,  we  obtain 


^3,/,=   a%=   a'lU=-U/.^      ,  (60) 


or 


>lh~    f^  Y  —  ^T 


^3/l/J  =  ^  ?,/>=-  i^»  ^c  ,  (61) 


> 
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which  is  Weingarten's  formula. 

From  Equations  (59)  and  (60) 


to.^  II  Y  -  O^tf  \\p, 


=      O 


(62) 


which  are  the  Mainardi-Codazzi  relations. 

The  covariant  base  vectors  of  the  space  (46)  can  now  be 
written  as  (16) 


9.=   '^- 


93=   ^3 


+  e  «/ 


b,°(- 


=     /^.  «>» 


and  the  metric  tensor  of  the  space  (46)  becomes 
where 


/>=  ^r'^^, 


r> 


From  Equations  (64)  and  (65) , 


/^ 


/"/» 


=  i.Vj' 


(63a) 


(63b) 


(64a) 
(64b) 
(64c) 


(65) 


(66) 


Raising  the  indices  Oi  in  Equation  (63a)  leads  to 


9  =  5    r«   ^/> 


(67) 


> 


or 
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f^=  g-'V- V^'=  ^r'^v^'   '  (68) 


where 


is  the  inverse  of  f^.      (16)  . 


2.4  The  Relationship  Between  Space  and  Surface  Quantities  in 
Normal  Coordinates 


Differentiating  Equations  (63)  with  respect  to   9    , 


3.,^=  A^Y'^-'r^" 


(70) 


From  Equation  (21) , 


r  -r   .  r-i 


(71) 


Substituting  Equation  (56)  into  (70),  and  combining  with  (71), 
Eliminating  9»  ^^^      3j  from  Equations  (63) ,  Equation  (72)  yields 

;'  (^:^^r-^l-r;;.t)%^(r^/i-r^)s=o  (73) 

where 

V      i  „»y   J        ^»>  1  (74) 

Fi  Equating  the  coefficients  of  O.^    and   a.j  to  zero,  respectively, 

r/^  =  (/*"/^(/^I(|/*  --  fJ^  )    ^  (75) 

r^^  =  y"x  r  (^    .  (76) 
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These  are  the  relations  between  space  and  surface  Christoffel 
symbols.  Equations  (75)  and  (76)  have  been  derived  by  Naghdi 
(16)  but  in  a  much  more  lengthy  and  indirect  manner. 
From  Equations  (23) , 


(77) 


Referring  1/  to  the  space  (46)  with  base  vectors  Uo.    .  a,, 

(78) 
Comparing  Equations  (77)  and  (78),  with  Equations  (63)  in  mind. 


-    V^  a  -^     V)  a.    -   V   ci^-^V  a.j 


■*"•*"  ~         .i,«r    fly 


>    ,»     , » > 


^3=  V*=  \/',  =  V 
Equations  (79)  are  the  relations  between  space  and  surface  vectors 

v/*3  V 

(16).   V   is  a  surface  invariant.  The  A^  and  its  inverse  act 
as  "shifters"  in  the  space  of  normal  coordinates. 
From  Equation  (28), 

which,  with  the  help  of  Equations  (65),  (76),  and  (79),  reduces  to 

y    »  ■  ' 

V«|^  =  /^«  (  V^/l^  -  b^V/")  .  ^81) 


mr 
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Similarly, 

,  V\=     C/--X  C/7J,-b:v'0         .  (82) 

For  the  other  components  of   Vi  Ij  and  V  Ij   ,  we  obtain,  in  a 
similar  manner,  the  following  expressions  (16) 

Equations  (81),  (82),  and  (83)  are  the  relations  between  the 

space  and  surface  derivatives  of  a  vector. 

-■  '  ~^.  ■■■-   ■ .  ,>• 

2.5  Modified  Hellinger-Reissner  Variational  Principle 
For  an  elastic  body  (18) 

where  ^   is  the  strain  energy  function,  measured  per  unit  volume 
of  the  undeforraed  body,  and  satisfying  the  following  relation 

Si=  s^^  SXj      '  ,v::  :  .     ^       >  (85) 

where  ^   stands  for  variation. 

The  modified  Hellinger-Reissner  variational  theorem  can  be 
stated  as  follows:   The  state  of  stress  and  displacement  which 
satisfied  the  differential  equations  of  motion  and  the  strain- 
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displacement  relations  (39)  in  the  interior  of  the  body,  and  the 
conditions  of  prescribed  stress  on  the  surface  of  the  body,  is 
determined  by  the  variational  equation  (14) 


i 


•-•'^  ..  .    ,  ,i  _       (86) 


=  O   , 


in  which  the  components  of  the  strain  tensor  VI   ,  stress  tensor 


'.S 

5"'''   ,  and  displacement  Vi      are  allowed  to  vary  independently. 
The  5.   and   ,t>   denote  respectively  the  density  and  the  com- 
ponents of  body  force  per  unit  mass  of  the  undeformed  body,  o  Ou 
the  components  of  the  acceleration  vector,   •'C   the  volume  of  the 
undeformed  body,  ,h  its  total  boundary,  where  only  the  stress, 

t  is  prescribed,  and  dZ     and  dS      represent  the  corresponding 

elements  of  volume  and  area,  respectively. 

By  Equations  (18),  (51),  and  (66)  the  element  of  volume  is 


) 


=  ^dA6e  ^^  ,  (87) 

For  the  face  boundary,  the  element  of  area  is 

^^=  /^  JA     .:  ,  (88) 

and,  for  the  edge  boundary  (16), 

■?  '  ■  ' 

J.   .  '  ,  ■   '  ■ 

•  n.td5  = /^.v.c//5  de^  ,  (89) 
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where  o^o.   are  the  components  of  the  unit  vector  normal  to  the 
intersection  of  the  undeformed  middle  surface  and  the  edge 
boundary,  and  d/    is  the  line  element  of  that  intersection. 
The  variational  Equation  (86)  can  be  written  as 

S  [  j(  s'' r,j  -  $)dr  -  /[  f(  v.,  .  v^u  .  v\  ys\^  ^ 

•r 

.5   ^ 


V  ■■ "  •<■• 


(90) 


By  using  (81),  (82),  and  (83)  we  can  formulate  Equation 

(90)  in  terms  of  the  "shifted"  components  of  displacement  V  i 

*i  ,-,-'.  ■■_:.■ 

and  V    as  follows.  .   . 

+  .njt   ^/^■^I'f  *(.>T,t'  -►on,  t   ^Vjjib 


(91) 


in^" 
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We  now  carry  out  the  indicated  variation  in  (91) ,  using 
Green's  transformation  and  Equation  (87),  and  also  dividing  the 
resulting  volume  integral  into  three  parts,  to  obtain. 

Si,  *n,-  a,  *  sj  =  0  .  (92) 


where 


,.».  -a 


+  ^  J./.^  (.b"-  .«-)]  ^  /^  +  k^^[  5'^ 

+  [/-  s'V  v,V + 1^  \/l) . ^  s'V  /-H  \^,% )  ] ,  J  -^ 

r 


(93a) 


+ 
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..  .  •  "  .  • 


+ 


-/*s'V'!3J^v;  .  „n,^„t"- 


(93b) 


>  U,=/   |s-.(||-M_,s,^^,,,,.  ^  ^,3^, 


r""  "  (93d) 

r 


Here  ,C   represents  the  intersection  of  the  undeformed  reference 
surface  the  edge  boundary,  and  ,fi       stands  for  the  total  area  of 
the  undeformed  reference  surface. 

Since  the  displacements,  strains,  and  stresses  are  regarded 
as  arbitrary  functions  and  allowed  to  vary  independently,  we  can 
set  Si,     ,  n,  ,  ?I,  ^  and  ^J   equal  to  zero,  respectively,  and 
thus  obtain  the  four  sets  of  fundamental  equations.   In  what  follows. 
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we  shall  apply  this  principle  to  the  three-layered  sandwich  shell 
in  order  to  obtain  the  corresponding  equations  of  motion,  strain- 
displacement  relations,  constitutive  equations,  and  boundary 
conditions. 


> 


.  '   ;  CHAPTER  III 

GENERAL  NONLINEAR  THEORY  OF  SANDWICH  SHELLS 

3.1  Change  of  Reference  Surface 

■''"'-  '     ■  -i 

/.  —     » 

Let  2  h  J  2 h   J  2 n   be  the  thickness  of  the  upper  facing, 

core,  and  lower  facing,  respectively.   From  now  on,  the  single 

prime,  bar,  and  double  prime  are  employed  to  designate  the  quantities 

referring  to  the  upper  facing,  core,  and  lower  facing,  respectively. 

We  choose  our  right-handed  general  convected  coordinates  &        as 

a  set  of  geodesic  normal  coordinate  system  situated  on  the  middle 

surface  of  the  core.   It  is  found  advantageous  to  choose  similar 

coordinate  systems,  "   and  O      ,  in  the  middle  surfaces  of  the 

upper  and  lower  facings,  respectively.   Then,  by  the  definition  of 

geodesic  normal  coordinates,  we  have  the  following  relations 

e=  0+z  =  0-f2^    a=,0='e    .  (94) 

where 

s 

Also  let   V^   ,   V^    ,   Vi   be  the  components  of  the  dis- 
placement vectors  in  the  upper  facing,  core,  and  lower  facing. 
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respectively.   By  using  Equations  (79)  the  "shifted"  components 

'.  .       —  < 

of   ''«:   ,   l^i    ,  and  \/c     take  the  following  form. 


v. 


v: 


V.  (=f^1       V 


V> 


ft 


V. 


=   v; 


(96) 


;s 


In  general   V^  ,  V;   ,  and  V^  are  functions  of  Q''   , 
and  could  be  expressed  as  an  infinite  power  series  in  terms  of  the 
thickness  coordinate.  But,  for  the  sake  of  simplicity,  we  shall 
assume  that  they  vary  linearly  as  follows. 


Vt=  U,  -V  ce^'z)Yi 


v:=   ii.+  e^t 
V;=  V,^    (e^-%)Y, 

Substituting  Equations  (92)  and  (95)  into  (97),  and  con- 
sidering the  continuity  conditions  of  the  displacements  at  the 
interfaces,  we  obtain 


(97a) 
(97b) 
(97c) 


U,  =  l^,  H-  ht,  ^\\ 


■u. 


.-a.-  i;f.  -^Ti 


(98a) 
(98b) 


Elimination  of  U.i  and  U^  from  Equations  (97),  by  using  Equations 
(98),  yields 
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v:=  ^i^  ^'% 


*  .«• 


V  =    V,-.   e'Y, 


(99a) 
(99b) 
(99c) 


where 


(100a) 
(100b) 
(100c) 


■s  ■_<  t 


Equations  (99),  imply  the  following  strain  distribution  for 
the  upper  facing 


-»  '-. 


•or} 
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X, 


(101a) 
(101b) 
(101c) 


Similar  distributions  for  the  core  and  lower  facing  corresponding 
to  Equations  (99b)  and  (99c)  can  be  obtained  by  replacing  the 
prime  by  a  bar  and  a  double  prime,  respectively. 

3.2  Equations  of  Motion 

We  now  write  the  volxome  integral  J  I,  ,  of  (93a)  for  each 
layer,  consider  Equations  (99)  -  (101),  and  integrate  with  respect 


M:-- 


> 
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to  U   through  the  thickness  of  the  undeformed  composite  shell, 
to  obtain 

S  I.  =      I'l,-^    il,  +  IX  ,  (102) 


where 


-f  v\  -  bV  ^)  +  '^^Vf\-bf  -  b"j^3)  -^ 

-'/^-JL;.-  'f'-'f^}rv,  .[C'-A^./.- 
^b^V,)  .  '^•'^^,.;,-b^1',)J,|.>['^'*'(  s;^ 

-^  '<  ->".v.)  -  '^^Vt\-b"-t-,^,)]^- 

^'C'    _'rn*]    §■%    j    dA  (103) 


'TT'S 
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^ 


Expressions  for  ? I,  ,  and  Si,  ,  can  be  obtained  by  replacing 
the  prime  by  a  bar  and  a  double  prime,  respectively. 

The  various  stress,  moment,  body  force,  and  acceleration 
resultants  appearing  in  (103)  are  defined  as  follows. 


'e-h 


ie' 


i-'h 


a- 


V'*= 


'.A 


'4 -'A 


't-'U 


n  =  /  yu  5  de 


"?-h 


C  =  M  +  C*  -  C, 


/u  's'/^cy.;  -  vi^-  b;^ .  eWi^  -b;  t,)J  + 


"■^ 


's^'t" 


e=z±'*. 


'pN 


/"'s'^r  ^.^  -  b^^«  +  eW^  -  ifV;)] 


(104) 
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->», 


■'•■-■   .  '  .  '■  -■'  "i  '    .'  .'-:■■ 

Substituting  Equations  (100)  into  (103),  and  setting  the 
volume  integral  SI,  equal  to  zero,  for  arbitrary  and  independent 
variations  of  the  twelve  unknowns  l^i    ,     ^fc    ,    ^%     ,  and  i'% 
we  obtain  the  following  twelve  equations  of  motion 


r  y 


(105a) 


> 


St: 


l%- 


32 

ujtj)  *  (W^'f:v^(t,.^.b;,'"(i)  *  (V* 

+  l;v^(V,.^*b*/^,)  .^.'^-t.M'f,* 

lU  *     '^''^  ^j  +  ^"Jf  -  bf  -^j)    ■*     - 

4[;^'»'.b(V*.n''^)J  Ctlj-bWj  -'n'»'t^ 

»- 

^C'm''-^;«'»*)('t%-b';t)   ^("m'*'- 

-hii'^'^ScYi.-brtj  .rt"-.^^'Y% 

-Vrjt,^^  ?'-  f'  =  0 

(105b) 

■[;F^-M;('n'"ivi-''jj (s'^uyk;;,^,)  -[l-^Ir-'•- 

-v-*)]  b^ .  r*'''-An''^v'')j  (tV-«>;;^.)  * 

4  h(w,'"L ;'„•"*}  ('^^11^-  b^'t,)  -  u ('w.-^ 

-];v'*)(TV-t;.Vo-bcY^'-vV)j,,,- 

-  ■f;;i-^ii(n-'LV'')j  (u,,^*b;u.)  .ll'''*.V(V^ 

-.*«•"*)]  ( t.^^  b^i)  ^  (i  ( - '^-  ^  "'"^ '  V ' 

-b;;t.)  - Tc-"^*'^ •""'*) ('^„;*-b;vo  - 

4[ro^^,)  -  Vci^t,)]  Ji 

■  !.:*<  t\  - 

_7;«j^^_  ?-'(5S:i\-blu,)    ^ 

*   C      -   nn       —     0 

(105c) 

^[^'^^Jcn'-n'")]  (U„^*bl,uf)  +£^-^rfV^ 

-^'k^*^;]  C  t,.^^  b^-V^,  )  ^  h  c^'''- *^ 'n'O  (  V.,^*  t^tj  - T  C'r^ 
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i 


(105d) 


4h*('n'«.V^J   C^lj-bltO   -li^'-lC'n^'' 
-'%)]    ^  'S     -    'd      =0      ,  (105e) 


34 


> 


► 


I 


i\ 


IX- 


^^x  ',|J  •'^'^  <"    ^^  -^>'  '     '  *^'"  *'^*" 


4  ];n''^)[i;-  ^x-f^('t\-v«f)  - 


(105f) 


(105g) 
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where 


-»-  cVl  -  W;  ( ,.  V3)  *  y  ^'d    =  o      ,  (i05h) 


n  '  =     n  '+   n   -I-  h 


^^=  '^^  ^^  V 


i 


I 


; 


C^=  M'^u'F^-r')  -~^('F>t) 


J 


J'=  _r'^'^-^'"■F^) 


> 


(106) 


ji^w  -;-^  :. 


> 


'      \      '  '  ' 
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3.3  Strain-displacement  Relations 

Writing  (93b)  for  each  layer,  considering  Equations  (99)  - 
(101),  and  integrating  across  the  thickness  of  the  composite 
shell,  we  obtain 

where 


V;  ,  ■.  V-  ■■   :         -tu^)  ^  (^C-b!^)C%-b^^>)  ^ 

■   -ir-bl(V,tf^-t,^^>)-fcii^V-i'f-to -^ 


:y-^^--v 
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> 


> 


► 


(108) 


Expressions  for  ^  I^  and  [  I^^   are  similar  to  Equation  (108)  . 

Setting  JIj^  equal  to  zero,  the  vanishing  of  the  coefficients 
of  the  arbitrary  and  independent  variations  of  the  stress  and  couple 
resultants  in  (108)  leads  to  the  strain-displacement  relations.   By 
using  (100),  we  obtain 

-  i>^  Sj  -  ui  c^j- '1-0]  f  ujj^  4  (;( tj,^  - 

4  X  ( -t,.,  -  ^,,- )  +  bl  uf  *  i;  bl  ctf  -  4pj  [  ^  Y  ^ 

'  _  ,  (109a) 

4    k     b^C+.-i.)]     j  ' 
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» 


> 


■ 


X/i=   2    I  -  t«  (  %li^  -  bj^  1-,  )    -  ip^  (  i'stt^  -  hg^  ^, )    -. 


Y    = 


.X, 


y„  =  ic^%*  t.r-ct,/] 


(109b) 


(109c) 


(109d) 
(109e) 
(109f) 


> 


> 
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The  strain-displacement  relations  for  the  core  and  lower  facing 
can  be  obtained  by  replacing  the  prime  by  a  bar  and  a  double 
prime,  and  letting  "  ""       and  h --"  ,  respectively. 

3.4  Constitutive  Equations  ' 

By  analogy  with  Equations  (45),  the  stress-strain  relations 
for  the  upper  facing  can  be  written  as 

V=  2  'C"*/*'  (  %»  -  Of^y'l  )  J  (110b) 


(110c) 


Similar  relations  hold  for  the  core  and  lower  facing,  respectively. 

Writing  the  volume  integral  d Ij  of  (93c)  for  each  layer, 
considering  (99)  -  (101),  and  integrating  across  the  thickness, 
we  obtain 


ti3=  i\^  n,^  a,      ,  ...        .  ail) 


where 


U,=    {[.„...,,$ -i|.)J^x,^[ 


in  - 


'/»   "  "/*" 


^  on 
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I 


and 

'   'l  •  ,  .  , 

is  the  strain  energy  function  per  unit  area  of  the  undeformed 
reference  surface   ©A 

The  strain  energy  function,   $   ,  for  the  elastic  anisotropic 
medium  subjected  to  a  prescribed  steady  temperature  field,   I    is 
now  assumed  in  the  following  form  (14) . 

^=  i  'S^nV.,  -  ^gl-)      .^    .  (113) 

*  •■  • 
Substituting  Equations  (110)  and  (101)  into  (113)  and  then  (112), 

^  and  setting  the  Equation  (111)  equal  to  zero,  the  following 

constitutive  equations  are  obtained 

+  ,0    A^   ».i        > 


(114a) 


(114b) 


iA 


/-•< 


^  =  2(  .8  '  .){;>  -^  .B   .^,  -  .T  )   , 

I 


(114c) 
(114d) 


?»   '   '     f  (114e) 


IP? 


> 


'n"=  :b''''x,  -r^'x^  ^'^'''y.r 


~)*r   '>*  ^ji    - 


where 


i  ♦  It 


'3-'*''=   ^ 'c-'* 'NeT  <^  e'   (n=o,,,z,3>; 


;6-'^»=  L  'c-'*"  Ca'/  d9^ 


(  W  =  O,  I  ,  2  J 


»-(. 

,'z^fc 


'B^'^lT/c.  f*^^0'r  d©'     (n«o,/,o 


.r"=/"l  t^'^'de' 


iA 


are  defined  as  isothermal  stiffnesses,  and 
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(114f) 


(115a) 


(115b) 


^-*'  (115c) 


(115d) 


^  -<j 


2-1 


«■-.--  '  -T 


> 


■'':\  '    -.  42 

as  thermal  stress  and  couple  resultants  per  unit  length  of  co- 
ordinate curves  on  the  undeformed  reference  surface. 

Similar  expressions  can  also  be  obtained  for  the  core  and 
lower  facing,  respectively.  ■  '.  ,  - 

■•',■>"  . 

3.5  Boundary  Conditions  '   i  '■        -  '.^ 

The  surface  integral  i  J  of  (93d)  can  be  split  into  two 
parts  ^ J,  ,  for  the  edge  boundary  surfaces,  and  Six  for  the 
upper  and  lower  boundary  surfaces,  so  that 

S  J  =  ij,  ^  ^Jx     ■■  ■  »  (117) 


where 


{. 


oC   -A 


(118a) 


(118b) 


» 
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In  a  similar  manner,  we  write  i  J,   of  (118a)  for  each  layer, 
substitute  Equation  (89)  into  (118a),  and  integrate  across  the 
thickness  of  the  composite  shell,  to  obtain 


II  ==   ^T.  *  il  ^  S'X 


(119) 


where 


a. 


'^  -y)-i:%)  ^Yr]]  rtr^^l'S*-   ov^f'-^'c^,;- 
n7i/.)    -'r-^^c'i.j^y/i'.)   *Y^'-'%^]]iV^- 

►  -fc>-)  -  V'c\j  -k>. )  * 

+  't''C'-1'0j  I  5^,1  cJxJ       .  (120) 


Here 


>:>- 


/iA 


■rHj^C'^^':'^'''!*']  ''^'     '        (ma) 


,~> 


,-'(■  =  /  A-„v.  i  '  ■!  .,  Me'  (121b) 


A -'A 


I  Expression  for  bj,  and  X  j,   can  be  obtained  by  replacing 

the  prime  by  a  bar  and  a  double  prime,  respectively. 

Substituting  Equations  (100)  into  (120),  setting  Equation 
(119)  equal  to  zero,  and  for  arbitrary  and  independent  variations 
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I 


» 


> 


of  the  displacement  components    6  M-i   i  '^i   S  t^   i  T^ 
we  obtain 

'•'■'■  • 

S  =  'S^  S    -'5      '       '"  .J-    ,  ^ 


(122a) 


(122b) 


(122c) 


w- 
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► 


» 


> 


=  .^,j 


('»n'^^.];'n'»^)[ij-  u//,.UtX-'tU- 


^^3  /^>  _     .A-J 


^»*: 


(122d) 


-  tj  u,  - 1 1)  ( t,-  ;t,) .  (V'- 1'-^^;  ['t'if  - 

4t-,1'.)    .t^.-t,)-/^ ''?'*(. -^^,)j       ,  (,22f) 

+  (^W]  J  (122g) 
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I 


► 


"e^  = 


"F.rr 


B.   Boundary  conditions  at  the  faces.   Evaluating  the 
integrand  of  Equation  (118b)  at  the  upper  and  lower  faces, 
respectively,  and  using  relations  (100),  we  obtain 


»        *  3 


/^X       rr.A 


'*  .r)-''(':-*p:)jsu,*-fc [<?>?:)- 


/^■^     »-,* 


/v/3         /"J 


(122h) 


OJA       *«A 


4(P.         ^         P.        )<Jt]     J         ^/» 


(123) 


where 


> 


/^'}'X 


'hA 


?:= 


'z-J. 


k/i. 


Setting  Equation  (123)  equal  to  zero, 


(124) 


O^t 


p:=  r: 


.  -p.  =?. 


(125) 


> 


^ 
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Here  the  signs  of  the  prescribed  stresses  are  considered  positive 
when  their  direction  coincides  with  the  positive  directions  of 
the  normal  coordinates. 

Up  to  now  all  equations  are  "exact"  in  the  sense  that  the 
complete  general  nonlinear  strain-displacement  relations  have 
been  used.  No  assumptions  about  the  state  of  deformation  have  been 
made  except,  of  course,  that  the  displacements  are  linear  functions 
of  0       .  The  "exact"  problem  of  composite  shells  is  governed 
by  the  equations  of  motion  (105),  the  stress-displacement  relations 
(109),  the  constitutive  equations  (114),  and  the  boundary  conditions 
(122)  and  (125).   These  equations  represent  twelve  equations  for 
twelve  quantities  ^i     ,    "^^      >  '^'^   »   '  i   ,  and,  thus  form  a 
complete  system  of  equations  for  the  problem.   In  this  form,  however, 
they  are,  for  most  practical  purposes,  too  complicated  and  approxi- 
mations must  be  introduced. 


» 


I 


I 


» 


CHAPTER  IV 


SPECIAL  APPROXIMATIONS 


4.1   General  Nonlinear  Membrane  Theory  of  Sandwich  Shells 

Equations  (105)  can  be  simplified  considerably  for  sandwich 
shells  with  very  thin  facings  of  identical  thickness,  "  .  For 
this  case 

^=  X=%  =  ^^=   o   .    5=5=0, 

and  the  last  four  equations  corresponding  to  the  variations  S  %  , 
oT}   ,   ST<»  ,  and   a  Ij  are  dropped.   The  results  follow. 

-^'?tj]t)^*f-f=o   .        ^      (126a) 


4« 
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■     •  -m'^'i.)-^rr}io   .f\f  =  0    ,  (126b) 


^ 


+  C     -  m^    =    0  (126c) 


+  (,!;>)+  CtX-"")(-'t,)  *  C-m'=  0    .        ,;       (i26d) 


'cfi,-;^' 


> 


> 
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These  are  six  equations  corresponding  to  the  six  displacement 
^.  functions  ^i     and   it  . 

The  related  boundary  conditions  can  be  derived  from  (122) 
and  (125)  as  follows 


and 


(127a) 


-^b'st)-^'*-^^t,  }  ,  (127b) 


(127c) 


(127d) 


(127e) 


A  further  simplification  of  (126)  and  (127)  is  possible  for 
the  case  when   i j  =  <?  ,  so  that  one  more  equation,  corresponding 
to  the  variation  %'^-^     ,  should  be  suppressed  in  (105).  The 
equations  of  motion  then  are 


r^  .,1^  r- 


► 


> 
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^  + 


I  -^f^-  f*=  0       ,  (128a) 

-^  1^  -  "P  =  °   .  (128b) 

-[  A•'^^;(w.■*'*- v")]  t^ .  [r'^-fi'c'^''^  v'')jx 


-^"(el^^'iu-C^,)  -ctiU-n  )t  - 


C^-  m^   =0 


(128c) 


and  the  boundary  conditions  are 

S'=  .>!»{^''C^^  -^^>1i^,)*C-'**-'i("'*'- 


(129a) 


fvr^i'^r- 


> 


> 
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:"rr  -'V^^^JlAt-^  ^^}    .  (129b) 

k  ■         '      :  .r('n'»'.V^Jt1,*[-fe'''fcr-'^''-*'J«'^tV}  ,      ^^29c) 

and 


/^/^      — t  O''-      — u 


P.  =   I*         .         -  Y-  =   y-         '  (i29d) 


t:  t    ■s'- 


where  the  definitions  of  P+  ,  r.   are  also  simplified  accordingly. 

4.2  Partially  Nonlinear  Theory  of  Sandwich  Shells 

The  modified  Hellinger-Reissner  variational  theorem  of 
Chapter  II  is  now  used  in  the  formulation  of  a  simplified  non- 
linear theory  of  sandwich  shells  on  the  basis  of  the  following 
partially  nonlinear  version  of  the  strain-displacement  relations  (39) 

^'j  =  T^^ilj*yj'^  +  V'I,V,|,)         .  (130) 

These  relations  had  been  used  by  Ebcioglu  in  (12)  and  (23)   for 
deriving  nonlinear  theories  of  plates  and  shells,  respectively. 
Using  (130),  the  variational  equation  (86)  becomes 

S  I  JC  5''  X-i  -  i  ]  dr  -ji  s''<Vi|j-  *  ^ju  *v'/.  v,|j  )dr  + 

+  j  ?.  M,'-  ,a^)v,dz    >  J."jt"^i  dS  {=   0        .  (131) 
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» 


» 


> 


With  the  help  of  Equations  (83)  -  (83),  we  can  shift  the  components 
of  displacement  into  the  reference  surface,  so  that    ■  f 

J  [  s'^  r,j  -  ij  dr  -  ]  i  S'''*C(/,)(  vr,i^-  b,^v;)  * 


>i 


z=   O 


(132) 


Executing  the  variation  and  using  Green's  transformation,  we  obtain 


where 


4  v;,»)]i)^  -^  [^  s"(  c  -  fe.1  v; )  -  /-  ^"(  I  ■* 


(133) 


(134a) 


54 


I 


> 


> 


r 


§      1  $ 


<i/»<je' 


<:'/)dfr 


(134b) 


(134c) 


•c  -I. 

-^v:.,)]  W;}  JA     ,  (134d) 

In  a  manner  similar  to  that  of  the  previous  chapter  we  obtain 
the  field  equations  of  a  partially  nonlinear  theory  of  composite 
shells  as  follows. 

The  equations  of  motion  are 

n;v.) .  f-^^r^.MH  *"'?"'''>  jbl- 


f  -  f  =  0 


(135a) 


:.'*        *  .  •     -(ff^^ 


» 


> 
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So:,:       j  n-^  a,.^-  b|.ii,)  ^  [--^^  UV'!.V'*;J(f„^- 


^■'^,-Vt,],,.    .   rn-^m-H/H. 


y^  -^ 


^i>'    -    i    =    O  ,  (135b) 


+  tl  V^  ■"^^""-'^'^^^'-''^j)  ^<^*-  "^^  =    0     ,        (135d) 
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I 


► 


I 


SI. 


St 


['% 


Ik.. 


(135e) 


'^■^      ',' 


■►('tll--n'')(<-tO-'3   -  a     =    0         ^  .     (135f) 

-   ^   -  3     .   ci     =     o        ^  (I35g) 

{(v.-^^  rn-^;[  ^3,^-];  (t,^- ^.^)  -  b^  «f  -  j;  (.^(-t, - 


where 


ft 


'c: 


-  '<¥ 


c^^^vn 


eO 


e'=i.'». 


(136a) 
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i 


I 


I 


^»='2ih   (136b) 


and  the  definition  for  f   and  'Kn**  are  based  on  Equation  (104). 
Similar  definitions  for  the  core  and  lower  facing  can  be  obtained 
by  replacing  the  prime  by  a  bar  and  a  double  prime,  respectively. 
The  strain-displacement  relations  for  the  upper  facing  are 

t 


(137a) 


(137b) 
(137c) 
(137d) 
(137e) 
(137f) 


Similar  relations  for  the  core  and  lower  facing  can  also  be 
obtained . 

The  general  constitutive  equations  for  r)       ,    >^       ,     "        , 
etc.  are  similar  to  (114),  keeping  in  mind  the  simpler  functional 


fT5»-r.  •   ^  ■  ■•■■  "'^\*  '  ■      .        ■         <r-  ■__'        ,_  -  .    -5  >   -  .•  ^.^^•■■f^'f 
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relationship,  expressed  by  (137),  between  the  components  of 
strain  and  displacement. 

Finally,  the  boundary  conditions  for  this  case  are  to  pre- 
scribe 


I  _  f.=  F.  .     -  F.=  P- 


and 


S  =  o>Jvf  "   -tjrn''  )       ,  ^38^) 

r=  o^Jn'*^(u„,.r,u.}  ^C^'*^-W'"'*-'"'*^)J(t,,f4t;t-)- 
> 

-hy%)  ,Y^r''^.*"'i''i,^^%}  ,         (138b) 

I  f  =  oy^[l^'^r^('n^-^'^)]-li'''^lC^'^'-'.^^ )][,]]       ,     (i38e) 


;       'e'=  .>i{('m'^'-j;v^)[u,,;-i;(t,.i-t„,)n;u. 
I  '  +KtWi-'i'«)]  -  (^'*'-r-''')(^,.f  ^ 


(138d) 
(138e) 


«^A 


*e   =  ,);,[(VMiV-)-(V^];v.'»*)fe»i] 


(138f) 
(138g) 


I 
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■^[,'shj   *Yi<*%)-~^Y('*%)}  ■  (138h) 

* 

4.3  Analogy  to  Donnell-Mushtarl-Vlasov  Approximation 

Simplification  of  the  above  equations  is  possible  under  the 
assumptions  discussed,  for  example,  by  Mushtari  and  Galimov  (19) 
and  Novozhilov  (20)  .   This  consists  in  neglecting  the  terms 
containing  the  product  of  ru  V^  in  the  expression  for  T-a  with 
the  following  results  for  the  strains 

jl  ^"S  =i  [/^-  V>  *  ^^''-^   *  ^•'  ^^*  *   ^'*-  ^'''  ^   '  (139b) 

The  variational  equation  (86),  compatible  with  the  nonlinear 

/■ 

St rain- displacement  relations  (139)  becomes 


or 
+  JC(.'>t'  ^  .n,t  )/^.V,  +(.n^t  -."^^t  ). 

.^     ^         1 

»  Vi  J  d  S   =   0         ,  (140) 


I 


» 


» 


where 


n3  =  r[s'-^-i(^-^)jsx,H^^^^  , 
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(141) 


(142a) 


/«r 


-o   ■^'((•i 


(142b) 
(142c) 


r'. 


«c  -A 


r  (142d) 


By  the  procedure  used  in  Chapter  III,  we  obtain  the 
corresponding  nonlinear  equations  of  motion  as  follows. 


Ju^  •• 


r  r 


(143a) 


'^^rf-^rr 


l> 


I 


n 


l/«* 
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(143c) 


-  pi")  (  •+ -ti)   *  C'  -  m'  =    o  ,  (143d) 

'J     'J     '1^ 
t('tX-V')(,.^o-  '5'-  'cJ'  =    0    ,  (U3f) 

ft;;     [(W.rn-').(V''>];"-''''H^J„«-VUl  - 

W     i"  ,     J  II. I 

-^-^-6=0  ,  (i^3g) 
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I 


and 


ClA3h) 


^  ^    *  d    =    o 

The  corresponding  boundary  conditions  are  to  prescribe 


(lAAb) 


l^^    .^^j[m'''*];('n''t"n/'i  -[^"/'-WV.VOjb:]    ,     (144c) 

J     '  (144f) 
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^ 


I 


> 


"e^ 


(144h) 


where 


,  t  ^ 


[  t;  J 


[t-'^'''~'^,,f''^%.f-'^'tf)'/i"'.'''l,)W  \ 


e'.ii'k 


and  \^      ,  Gt  are  defined  in  Equations  (136). 

4.4  Partially  Nonlinear  Membrane  Theory  of  Sandwich  Shells 

When  the  facings  are  so  thin  that  they  are  effectively 
membranes,  we  set 


and 


5=5=0 


in  Equations  (143)  and  (144).  Accordingly,  we  may  suppress  the 

A   '  '  •  f 

equations  corresponding  to  the  variations  btj  ,  ft,   ,  ^  if  , 
and   ii,  .  The  resulting  equations  can  be  simplified  considerably 
for  sandwich  shells.  The  equations  of  motion  are 


(145a) 


(145b) 
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I 


-  '^   ^  '  -  "  =  ^      >  (145c) 


(145d) 


The  edge  boundary  conditions  are 

4.5  Sandwich  Shells  vith  a  Weak  Core 

Assuming  that  the  components  of  stress   S    in  the  core 

I  a^e  of  negligible  importance,  they  may  be  set  equal  to  zero,  and 

the  components  of  transverse  shear  and  normal  stress,   S*   and 
—  >j 
5     ,  only  are  retained.   The  equations  of  motion  (145)  and 

boundary  conditions  (146)  then  become 


^/,'^:^^', 
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I 


I 


+  (rt- ?;";(«*  t.)  -  c*  -  '^^  =  0 


(147d) 


and 


(148a) 


+  ^''c  '  -^-t,)  j    J 


(148b) 
4  Pc  i^^,)j    ,  . 


(148d) 


respectively,  where  the  definitions  of   Kl"   ,   *^"   ,   f    , 
and   T   are  also  simplified  accordingly. 

.       r 

4,6  Linearization  of  the  General  Equations  of  Motion 

Dropping  all  nonlinear  terms  in. the  equations  of  motion 


p-'?.' 


^ 


I 
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(149a) 
(149b) 


(105)  and  in  the  boundary  conditions  (122),  we  obtain  the 
following  equations  of  motion 

fa,:    f «.  *  (  n'*"-  ^"^  ^^  )  ^  -  p'  -  f  *  =  0    , 


+  C  -  m  =  0 


(149d) 


/-r  /' 


-  ^     ^     3      *    d       =      0  ,  (149g) 


and  boundary  conditions 

'  ■   (150a) 


S'=  .>iCn^-r.''^t;) 


ifimsm'' 


» 


► 
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r-  .v./ 


>    '  >  (150b) 


(150d) 
(150e) 


"e'=    .>^[(^-''^V'0-(V'-i;V^)i>;]      ,  (i50g) 


CHAPTER  V 
COMPARISONS  AM)  CONCLUSION 

The  nonlinear  fundamental  equations  of  a  theory  of  sandwich 
shells  in  terms  of  the  undeformed  state  have  been  obtained  from 
three-dimensional  continuum  mechanics  with  the  help  of  the  modified 
Hellinger-Reissner  variational  theorem.   These  include  strain- 
displacement  relations,  equations  of  motion,  boundary  conditions, 
and  constitutive  equations  for  elastic,  anisotropic,  composite 
shells  subjected  to  large  displacement  gradients  under  the  influence 
of  mechanical  and  thermal  loads.   Several  approximate  systems  of 
equations  which  may  be  suitable  for  application  to  cases  in  which 
the  displacements  and  rotations  are  restricted  in  magnitude  have 
also  been  obtained.  . 

For  a  sandwich  plate,   M*   is  simply  the  Kronecker  symbol, 
since  then  o»\  ~  0   .   Equations  (105)  become 

+  Y  't,  ^  ^^  t,  *  Y  't)  ]■  1^  *  f'-  f '  =  0     ,  (151b) 

68 
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*      -i. 


<;  ^,  :  j  [^'^^  I  Cn-'*-  "n"^)]  (  h^  *  a^^)    -  [  :j^'^  h^'.-^ 


.Vo-o  Vi^  .];(fY'-ry)]l^-(^L- 
^  -  h*\  i  +  't, )   +  C     -  m    =    0        , 


(15Ld) 


1-3  -  d    =       o 


(151e) 
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l\ 


1% 


4cV^^iW'*)Vl^*"^''1•'V|l«  -^   . 

-    <]     -  d  =  0 


(151g) 


(151h) 


> 


It  is  easy  to  see  that  our  results  can  be  reduced  to  those 
obtained  by  Ebcioglu  (11) ,  and  shown  by  Ebcioglu  to  contain  the 
earlier  works  of  Eringen  (21)  and  Yu  (22) . 

For  homogeneous  shells,  we  simply  drop  the  upper  and  lower 
facings  of  the  composite  shell  and  retain  only  the  core  layer. 
Then,  all  terms  including  quantities  denoted  by  single  and  double 
primes  and  the  last  four  equations  of  Equations  (105)  should  be 
suppressed.  We  have 


i     T 


7* 


(152a) 


"fjgp-- 
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+  p  -  f  =  0   , 

^'tj :       j  m-^  u»,^-  bj.^^)  ^  r^  t,,^-t^'t.)  jiu  ^ 


(152b) 


(152c) 


(152d) 


I) 


This  is  in  agreement  with  the  results  of  a  nonlinear  theory  of 
elastic  shells  obtained  by  Habip  (_13)  ,  (14).   In  the  absence  of 
curvature  effects,  these  also  reduce  to  those  of  Habip  (13)  for 
plates  in  the  reference  state. 

In  the  partially  nonlinear  case,  our  Equations  (135)  agree 
with  the  results  of  Sanders  (24),  given  for  the  small  strain 
approximation,  if  we  drop  the  effects  of  the  facings  and  adopt 
Kirchhoff's  hypothesis  for  the  resulting  homogeneous  shell. 
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» 


> 


By  setting   'j^=  J,  =  ^5  =0    ^  and  noting  that 


Equations  (143)  lead  to 


n'-'L      -  fbS  P^-  f^  =  0   , 


i   _> 


(153a) 
(153b) 


7   ,  '/1«     */!•'  »  I  i 


\%:  [  ^"^  I  ro-^-  "n"')]/!.  -  W  'f  -  r  )!.:  - 


r%: 


w, 


*    ^    -  d    =     o 


(153c) 


(153d) 


(153e) 


These  are  identical  with  Ebcioglu's  equations  (12). 

For  comparison  with  the  well-known  works  of  Reissner  (6), 
(Z)>  we  transform  Equations  (147)  into  their  conventional  form, 
in  terms  of  physical  components,  as  shown  in  the  Appendix,   In 
order  to  show  that  these  previous  results  are  contained  in  (147) 
we  use  the  following  notation. 


9,  =  X      ,      e^  =  y     ,    /«:,=, 


'flx,  = 


,    R,=  R.=  «' 


N..=  M..=  Nx;   N.,=  N..=  Mx.j  ^   K=N„=Mj  , 


a,=  Q, 


a.=  ds       .    M„-  /iM=M« 


ArN,ra,    .    M„=N.,=  M5,,    P.=  P^=c,  =  c;,=  o   , 
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\ 


and  also  omit  the  effects  of  thermal  gradients  and  acceleration 
resultants.  Equations  (A8)  bectxne 


Qx 

2^    , 

1^=      0 

rix 

I'd 

^^^    . 
'i'3 

9  X                 3 

-^'^s    9- 

-M,    -^e      ^p     _^ 

^(^' 

+  4  c  N. . 

..^/«    .^^/..  ^'e      , 

2  AJ.. 


?  v^ 

i^u 


3   -^x  73 


M 


r^ 


-;  -^ 


-  crj„  (  I  +  e)  +5  =0 


(154a) 
(154b) 
(154c) 
(154d) 


(154e) 


(154f) 


\ 


These  are  the  quations  given  by  Reissner  (^) ,  (2). 

Completely  linearizing  Equations  (A8)  and  neglecting  the 
acceleration  resultants,  we  obtain 


9  0. 


")  Q. 


ff^N„./^(-f  *P.}=  0   . 


(f^  N,J 


a»N,x) 


9/«T, 


1 

3d,* ^^'    30/  ^     .    ^  ^^ 


A/»,  - 


^e, 


«, 


(155a) 
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A  . 

.  _^«  )  -^  Jo-  P  =   0    ;  (155c) 


|-       '     -^,(/^^'^..)--|^//^^^.)-  feT^ 


y^^ 


M^.^J^{C,-Ci,)=     0  ,  (I55d) 


(155e) 


^  These  are  the  equations  obtained  by  Reissner  (5)  for  the  small 

bending  and  stretching  of  sandwich-type  shells. 

'.  '  The  above  comparisons  reveal  that  the  modified  Hellinger- 

Reissner  theorem  is  particularly  suitable  for  the  consistent 
derivation  of  intermediate  theories.   However,  during  the  course 
of  this  derivation,  we  found  that  the  order  of  the  variation  and 
"shifting"  is  important.   If  we  execute  the  indicated  variation 
in  Equation  (86)  before  shifting  the  components  of  displacement, 
i.e.,  before  we  employ  the  normal  coordinates  in  the  three- 
dimensional  body,  some  inconsistency  will  appear  between  strain- 
displacement  relations  and  the  equations  of  motion,  forcing  us  to 
drop  certain  terms  in  the  strain-displacement  relations  in  order 
to  make  them  compatible  with  the  equations  of  motion  (14) .   This 
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arises  only  in  partially  nonlinear  theories.  For  the  general 
non-linear  case,  the  order  of  variation  and  "shifting"  is  im- 
material. 


» 
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APPENDIX 

For  the  purpose  of  comparison,  we  now  proceed  to  record 
the  physical  components  of  the  equations  of  motion  (147)  for  the 
nonlinear  membrane  theory  of  sandwich  shells  with  a  weak  core  and 
identical  facings,  i.e..   In  =  f)  =  t    .   The  material  of  the 
composite  shell  is  isotropic,  and  the  coordinates  are  in  lines  of 
curvature.   The  body  forces  are  neglected. 

Following  Ebcioglu's  procedure  (12) ,  assuming  the  membrane 
and  transverse  shear  stresses  to  be  uniform  across  the  thickness 
of  the  facings  and  the  core,  respectively,  and  with  the  transverse 
normal  stress  given  as 

S"=  ^-  ^  ^-  .    J^(%.-^..)         ,  (AD 

we  obtain  the  stress  and  couple  resultants  in  terms  of  physical 


components  as  follows 


h  -t 


'rio.)  =  I    '3c,.,  c*  ^J  ^  ^J  =    N.,  (  .  *  -^j     > 
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-''-li 
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(A3) 


=  2I  cj;,n 


where 

and  ^jm   represents  the  value  of  the  effective  transverse 
normal  stress. 

From  Equations  (A2) ,  we  can  deduce  the  following  definitions 

n(..)  =  'n(,o  ^    "<•-)  = 


Hoo  =  '•^(.u  *  "'^(")  =  '^.-  *  ''^"  =  '^" 


n,*l,  = 


=  C I  +  — r— ;  Nil  +  ( I  — 5 — ;  Nxi 


n<,»,=  'n<,»,  *  n,„)  =  N.t  *  '^^.t  -  '^1 


'■•i^, 


"'TT;; 
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•^(..) 


=  (I  -^4^)  K  -^('-^)M.,  =  N.,     , 


Ra 


R» 
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:,!)=   'Hc*,  *"»^(.x)    =   N,a-*-*N,i  =    M,t  =  N>,     . 


n(x,)=  ''^<»'J  ■"  '^<»"  = 


=  (,.1^)N,,  .  (,-i^)K=N.  ' 


l.-t 
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In  the  same  way,  the  following  effective  external  force  and 
moment  intensities  are  obtained 

C„,=  I[(-iif)(-^)L,  - 

-where   f+ii,)  and   P-a)   are  the  prescribed  loading  at  the  upper 
and  lower  facings,  respectively,  and 

If  we  let 

where  6  represents  the  effective  transverse  normal  strain  for 
the  composite  shell,  the  equations  of  motion  (147)  then  become 


n-^ 
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where   J^^  ,  T3  ,  ifCw  ,  <rUj  represent  the  physical  components 
of  the  inertial  forces  and  moment  resultants. 
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